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Запропоновано аналітичні моделі гравітаційного поля та геологічного середовища для горизонта- 
льно-шаруватого середовища з кількома густинними межами поділу шарів у нижньому півпросторі. 
Перша модель являє собою рівняння сили тяжіння, отримане з розгляду нормального потенціалу сили 
тяжіння у локальному околі поточної точки вимірювань у двовимірному зображенні. Для моделі сере- 
довища, на додаток до відомих в теорії інтерпретації потенціальних полів основних класів контактних 
поверхонь (Нумерова, Остромогильського 1 Страхова), введено новий клас контактних поверхонь -- 
клас Чорного. Наведено зразки теорем розділення полів для випадку кількох однозв'язних об'ємів 1 кі- 
лькох шарів, що не перетинаються (два найпопулярніші початкові наближення середовища) у дослід- 
жуваній моделі. Описано новий чисельний підхід до визначення початкового наближення густинного 
контакту і його асимптот у багатошаровому геологічному середовищі за допомогою кількох 'контакт- 
них" алгоритмів. 
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Основное направление научного прогресса -- 
зто создание общих математических моделей в гео- 
логий, геофизике и в дальнейшем - создание геоло- 
го-математических теорий... 

Из докл. "Матем. моделирование в петрофи- 
зике при комплексной интерпретации данньх ГИС, 
керна и опробования", ПАНГЕЯ. уммм.ММЕІЇ. пагоаї ги 


Введениєе. В последнее время большое внимание уделяєтся вопросам переинтерпретациий 
картографического материала, полученного по данньт сьемок 30-летней давности, ввиду 
резкого ограничения обьемов геологоразведки. Однако вследствиєе повьшенньх требований к 
точности и зффективности интерпретации такая переобработка больших обьемов геофизи- 
ческих данньх требуєт новьїх подходов к численному обоснованию методов, используємьх 
при подготовке и обработке соответствующего картографического материала. Вполне логич- 
но обьединить развитьєе математические конструкции с современньтм интерактивньм програ- 
ммньм обеспечением, дабьт комплексно решить проблемьт оцифровьгвания, визуализации и 
истолкования геофизических данньх. Проблема истолкования данньх, наиболее важная с то- 
чки зрения приложений (впрочем, тесно переплетенная на уровне первичной обработки с пер- 
вьіми двумя), решаєтся в последнее время в русле перестройки аналитической и вьштчислитель- 
ной базь  теорий интерпретаций потенциальньх полей, обеспечивающей становлениє новой 
(третьей) парадигмьт в теории и практике интерпретации. 

Современнье методь обработки геофизических данньх развиваются для создания цифро- 
вьїтх аналитических моделей поля и геологической средьт и вьісокоточньх алгоритмов, ориен- 
тированньх на применение в пакетах ГИС, особенно при отображений входной информации 
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и результатов есе трансформации с возможностью сопоставления с априорной информацией, 
также представленной В графическом виде (так назьваєм. визуальная интерпретация, рис. 1). 
Фо ї 
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ретро 5 
Рив 1, Карий Равнтаднониніх аномалий Украинь! в среде свободно распространяємой ГИС Опапішт О15. 
Собственно, ,послойная" обработка представляет новьшй зтап качественной интерпрета- 

ции, в котором возможности представления информации в виде злектронного журнала изме- 

рений |Якимчик, 2010) можно использовать для предобработки информации в среде специа- 
лизированньх пакетов общего и специального назначений. Примером такой предобработки 
могут бьть процедурьг: вьіделения регионального фона, например, путем минимизации невяз- 
ки между набором данньїхх и некоторьтм полиномом в среде Мацаб; предварительная фильт- 
рация измерений в среде Ехсеї; нахождениє особьгїх точек поля методом обратной интерполя- 
ции полинома Лагранжа в Маріє; редукции поля с помощью алгоритмов В.Н. Страхова, реа- 
лизованньве в Руїоп, и т.п. Одним из наиболее удачньх примеров цифровой качественной 
интерпретации можно считать обработку значений поля сильгт тяжести в связке ПО Маріпіо - 

СогеЇДгам, предложенньгй в работе |Якимчик, 2010|, что, отнюдь, не исключаєт их заменьі в 

обозримом будущем на ореп-50игсе аналоги. 

Коснемся, однако, иного направления в рамках детерминированного подхода к интерпре- 
тации данньїх потенциальньх полей -- разработки новьїх аналитических конструкций геофи- 
зического поля и геологической средь. Относительно математических конструкций поля от- 
сьлаєм заинтересованньх читателей к работе |Якимчик, 2009|, а среди моделей средьт авто- 
рам наиболее близки плотностньєе модели типа контакта, одну из которьжх предлагаєм внима- 
нию читателей. 

При разработке таковьх конструкций ради повьшения точности интерпретации данньгх 
геофизики обязательно следуєт учитьвать условия адекватности постановок современньтм 
требованиям к получению исходньїх геофизических данньжх и их обработке. В целом, суть 
зтих требований сводится к таким пунктам: 

е вьіокая точность и технологичность алгоритмов и программ, их переносимость; 

е- адаптация моделей поля к заданию исходньтх измерений на коротких профилях; 

е адаптация моделей геологической средьт к ее сложному (нелинейному неоднородному) 
состоянию. 

Определение контакта в известньїх моделях средь Страхова, Нумерова, Маловичко |Чор- 
ний, Дубовенко, 2002| традиционно замьткаєтся на вьтчислений разности ундуляций некой те- 
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оретической модели средьг от некоторого опорного ее злемента, от которого зти ундуляции 
не слишком сильно отклоняются. Подвох в том, что при зтом необходимо заранее знать неко- 
торьтге злементьт геометрии (глубину, форму, ориентацию) данного опорного злемента (моде- 
льного тела либо контакта), как начального приближения итераций. От вьшбора удачного на- 
чального приближения зависит, в целом, зффективность любого из существующих итераци- 
онньх алгоритмов |Старостенко, 197481. 

Предлагаєм вьведенньге для горизонтально-слоистой средьт с несколькими плотностньшїми 
границами раздела слоев аналитические модели поля и геологической средьі, в которьгх зти 
опорнье злементьт рассчитьваются в процессе итераций. Модель поля представляєт собой 
уравнение силь тяжести, полученное из рассмотрения нормального потенциала силь тяжести 
в локальной окрестности текущей точки измерений в двумерном представлений. Размерьт са- 
мой локальной окрестности существенно зависят от размеров самой тяготеющей области и 
точности аппроксимации самой нормальной составляющей потенциала | Дубовенко, 2011. 


Основное уравнение поля. Рассмотрим в трехмерном неограниченном евклидовом прост- 


ранстве Є (3) 


две области С, и С,, вложенньюе друг в друга так, что расстояние между их гра- 
ницами СС; и СС, удовлетворяєет неравенству азс. дС, ) 2 0 при С, СО, «6;. Рассмотрим 
в области С, начало координат и точки области С" «С, 26», С, з С 0С», будем обозначать 
через і з (В, В ба а точки замкнутой неограниченной области С" - ЕС, дС'з ОС, - 


через х - (ху дода). 


Пусть область С, заполнена массами с плотностью 5, (2), а область С, з СС» - масса- 
щ- з 
ми с плотностью 5 (5) (рис. 2). Обозначим через п(х) единичную нормаль к поверхности СС", 


внешнюю по отношению к С" ; тогда нормальная составляющая силь тяжести равна 


чув чізьоазсзіьцоо! 


Рис. 2. Модель контактной средь! в классе Черного. 
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| дУ(х). ду "(з 5) 
цх)- стою ) бекон а) 


где 5(5) з 5(2)- 6, (2), ЕєС», 52) «0, 6 є С, / - гравитационная константа. Вообще гово- 


ря, силу тяжести описьваєт функция 


их) є гад У(х) є / | 5(5)жгад "(х 5)а5, (2) 


поскольку значения силь тяжести - зто значения модуля градиента потенциала, а не составля- 
ющие притяжения, в том числе и (1) |ДДубовенко, 2009). 
Упростим вьгражение (1), рассмотрев локальную окрестность С точки х при большом ее 





зувеличении". Тогда с известной долей приближения |Дубовенко, 2011| участок ДС - ОпдС" 
границь ОС" будем считать ,неограниченной" плоскостью, т.е. 
ОС з їх зо надо о) з соп5і, (5 зе Я 
где 5, - С 00 0С" . Расположив в зтой плоскости прямоугольную систему координат так, чтобьт 
ось ох, бьша перпендикулярна к плоскости, получим 
п(х)є є, х 2 (ху Хода), 52 (б»бо»ба), 0" -О0СЯ, 
о еогсо,дасаОгодс. 


Пусть в окрестности 0 функции 5. (5) з соп5ї, бСеС,С0, 5, (2) зсоп5іі, СеО СС, 
и 6-8, -б, а С, - собственноє подмножество С. Фактические размерьт области О 


зависят, прежде всего, от размеров изучаємой тяготеющей области С, и точности аппрок- 
симации вьражения (1) | Дубовенко, 2011), т.е. в нашем случає 








їй Ї па, б Ї ЗБОРУ «є. 3) 


Если є» 0 достаточно мало, то (2 вьїюбрана удачно. Границу СС, - СОС», вьбираєм 
достаточно гладкой, 
ОС, з і да РАРЯЮРНОЇ (ху )є 25, і,.,К - 1,2,3; 105, « дСьух, є си, 
При таких допущениях о среде получаєм потенциал притяжения 


«аа «Бабзав, б 
,2з дхз о дз У -х,) й 


где Ко /О. 

Если массь, сосредоточеннье в 0 , генерируют аномалий и(х,), вьтянутьтве, положим, 
вдоль оси Ох» , то их лучше всего измерять по профилям, параллельньм 06,. Но формула (4) 
малопригодна для обработки таких измерений. Проинтегрируєм вьтражениє (4) по 5., т.е. 


вьчислим интеграл 


уау 
2 2 3/2 4 
20 (т ву ) 
2 2 
зу", 45, з у. Делая замену переменной ин- 
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5 ЗЗАРУКЬ ЗОК З 
,Фаагіє ,т -ТтієФ- 5 -» ПоЛУЧИМ 
Т со5'Ф 





гдеиебуохьУЄі» хуЄЕо ти 


таФ 


тегрирования у - тіЄф, дфу-  , 
05 Ф 
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Я л4ф/со8'Ф УТ у .. Р 
Ци,м)з уг| оо з ; Гсовдар о у віпф з 
2 /со8 р) Рой т Р 

ум З 7 93 7 


- "зум" чад ву" е2 -х чи ру о -х) 


Если о -ЕФ зззш'зум , иуу'Є 7, то рассмотрим предел Ци, ут) при 20) -» 00, 29 -5-00 














(что имеет смьісл при шизучений  пространства вблизи поверхности ОС, причем 
дізк(5 ух )«« Чіатд ). Так как 
і У з У 
Пт, з т 


І, 
, 122 2 | 212 
уоайроьу о ПОН 


І) 





у 
то Ци, м) з 7 - при 50) -»00, ЕФ -» 0, и, следовательно, 
и" У 


537 
их, х,)є 2К "2 оо Гару 
, Ге. а -- (8, .І яю 
где 5 - проекция С, на плоскость ох,Х,, К - /о. Положим х, -0, 8, - 25, 7 а, В 


тогда 


о Ь ор й 
Ї 253 458. -кПо ( 1 ж) бі 53 Ме , (5) 


5 (8, -х) фе, а (я "М у ЕЕ) 


где Зм б ее) и а а ЗЕ), а «56, Р - кривьге, ограничивающие снизу и сверху тяготе- 


и(с0) - ка 


ющую двумерную область 5. 
Модели контактньтх поверхностей. Рассмотрим несколько популярньх классов контакт- 
ньх поверхностей, для которьгїх гарантирована единственность решения обратньтх задач. 


Класс АХ. Остромогильского Озі"а,Б). Доказано |Остромогильский, 1970), что в зтом 
классе однозначно определяются переменная плотность З 2.) и граница раздела сред х, (2 ) 


по значениям внешнего поля. 
Граница оС раздела сред с переменньми плотностями задаєтся уравнениєм 


М -пябо(), 2є Іа, 6| 
хе) є - ні кеьі (6) 
причем дх; (а) збх, (ь) «0,х, (є, ) - кусочно-непрерьвная функция и 
бо (о)2 с (х-а) ,асх «ачо,а»0 ай 





жк) 4-х) ,ь- бах сь а» 

Подклассь Оці) и Очі класса Озі"а, Б) можно считать моделями антиклиналей и 
синклиналей соответственно. Все греческие буквьт - прямкіє! 

Класс В.Н. Страхова 5Ї (а, Б) контактов рассмотрен неоднократно |Страхов, 1970; 19- 


74а, 19746; 1975а, 19756; 1976). Граница ОС раздела сред с различньми плотностями задаєтся 
уравнением 


| -й, бл (Є,), 2, Є Га, Б) 
чен Е зі 


При зтом бо (а)-0, ть 20, жа 0 - кусочно-непрерьтна и удовлетворяєт на концах 


(5) 
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интервала |а, Б) условиям 
бо) за(5 -а)",асбхачб,а»д 
|) 26(ь-а), ь-б«2хб,6»0 


принимая на |а, Р| конечное (счетное) число зкстремальньх значений, так что множество 


ем щої непусто. Здесь 


Їа,в) 





(9) 


25 3 сі) з вирехіт (ві), - 25 Є ху з оз) іпежиг 25) ог, є, (10) 


где с, с, Є|а, Р|, а точки с,, с, могут бьшть либо внутренними точками отрезка |а, Р|, либо 


совпадать с его граничньми точками. В последнем случає при условий, что й, « В, , ср ас, Ь 


Если (а, Р) - бесконечньй промежуток, т.е., а -20, Б -- чо, уравнение границьт (8) и 
условия (9) принимают вид 


сес()е-по ба), 8 є(-бу, бо), 
«2)-в|чаЕЇ", Кіа, (1) 
«су-ь|«ВИЇ, р»о, 
где й,П20,а г, ех(су) ки ес І са Че.)- іні ем (8 у при су с, ЄЇ 4, 5, 
причем для определенности положим - 2, «Й, - 2, - Л,, 0 «П, «Й. 


Для контактной границь раздела слоев с постоянной плотностью установлена единствен- 
ность |Чередниченко, 1978а; 19786). 


Пусть С" - горизонтальная полоса, ограниченная сверху и снизу прямьми 2-й, 2-Й», 
0 «Й «Й, «чоо, и пусть С  разделена на два слоя с плотностями д, и б, бесконечно глад- 


кой кривой ОС", звездной относительно бесконечно-удаленной точки, которая имеет асимпто- 
той прямую с «ЛП, п, «Пп «Й,, причем 


с(2)-н«кЯ, а»о, |Дь»о, «Дедс (12) 
Тогда, если значенис й и величина скачка плотности Лб ша б, -д, ИЗВеСТНЬ І, ТО граница Ге 


раздела слоев по внешнему полю восстанавливаєтся однозначно. 
Доказано |Страхов, 19746), что решение обратной задачи о контактной поверхности класса 
(12) единственно, если известньт скачок плотности Дб и конечная дуга кривой 2(4) из (12). 


Введем специальное обозначение для класса (12): СВ оо, со). Класс (11) будем обозна- 
чать 510/(-с0, со), Очевидно 


Озі"З(а, Б) с 5179 (а, Б), СА? (-оо, со) є 519723 (-со, со) 5109 (а, Б) с 5109 (-с0, со). (13) 
Проблема разделения полей. В практике гравиразведки найбольшесє значениє имеют зада- 
чи, в которьжх внешнее поле создают сложньє6е распределения массе, например, п попарно непе- 


ресекающиеся обьемь С,, К - І,2,...,п, или п попарно непересекающиеся границь раздела 


слоев различной плотности. В задачах со сложньїмми распределениями масс единственность ре- 
шения обратной задачи верна лишь в довольно узких классах, т.е., при наличий большой апри- 
орной информации об источниках. Теоремь единственности обратной задачи гравиразведки 
ДЛЯ СЛОЖНЬХ распределений масс основань на теоремах разделения полей, которьтре позволяют 
свести решениє обратной задачи по суммарному внешнему полю п обьектов (рудньтх тел, гра- 
ниц раздела слоев) к решению обратной задачи для отдельньваох обьектов по полям от зтих обье- 
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КтоОв. 
Пп ример 1. Пусть внешнее поле создаєтся массами с носителем в виде обьединения ико- 


нечньїх попарно непересекающихся односвязньх обьемов С, , К - 1, 2,..., п, с границами дС,, 
не имеющими попарно общих точек. Плотности масс в каждом из обьемов С, предполагаются 


ограниченньми. 
Теорема разделения полей. Если известнь внешние относительно друг друга гладкие замк- 


нутье конечнье поверхности 72.,, К - 1,2,..., п, строго внутри каждой из которьїх 5, нахо- 
дится одна и только одна компонента С, носителя масс, то по заданному суммарному внеш- 


нему полю (вне обьединения обьемов С, ) поля компонент С, определяются однозначно. 
Из зтой теоремь легко получить теоремьт для сложньх распределений масс, аналогичньве 
теоремам Новикова и Раппопорта-Сретенского, - для обьемов С, , звездньїх относительно за- 


данной точки |Чорна, 2003), либо обладающих средними плотностями. 
Пример2. Пусть в п-слоистой двухмерной среде каждьй из слоев С, имеет постоян- 


ную плотность д, . Границьт раздела слоев 2С, - бесконечнье кривье, имеющие горизонталь- 
нье асимптоть 2-й сей, (А, 20, й, 2 0) слева и справа (рис. 3). 

Теорема разделения полей. Если каждая из границ раздела ДС, содержит горизонталь- 

ньй отрезок ОС, (Є « соп5і, а, 2 5 «Б, ), известнь проекция некоторой его части на ось Ох 

х 2 х 2 В, и прямье у, разделяющие границьт слоев, то по заданному внешнему (при с «0) 


суммарному полю ДЯ с точностью до постоянньх определимьт поля Де , отдельньх границ. 





2-Ї 4С, С, зб 
По зчаузчззвіьзчо: сзозрз з оозьзьоі2Уі 


Рис. 3.К разделению полей от нескольких слоистьтх источников. 


Следствиє |Гласко, Страхов, 1977|. Если при зтом известньт положения асимптот 2 - й. и 


2-Й; и скачки плотностей До ,, границьт оС, (Кк -1,2,..., п) раздела слоев восстанавлива- 
ются по заданному суммарному полю Дє однозначно. 
Система непересекающихся контактньх границ и плотностей с(С) восстанавливаєтся 


единственньм образом І|Филатов, 1974). 
Рассмотрим теперь решения раки интегрального уравнения 


т ху а) М 
, 5 
Й! -хузаД не -ю 


где 2-2(с), ас хр -"верхняя", 2-2, и ас аб -"нижняя" границь плоской обла- 





2(х,2) с 
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сти С, . Положим в (5) 2(х) - К 'М(ю), 2(20) 2 Н есопі, 2.(2) - «(2), где 2(2) принадле- 

жит одному из множеств (например, (13), либо другому) допустимьх функций (линеалу), тогда 
получим нелинейное интегральног оо по отношению к искомой функции 2 - 2(2): 

2 
або Ї 7 У ЯН 

солі нн 

Его можно рассматривать как зірвано на некотором линеале. Найдем его приращениє 

Ле(х, 2) - є (Х, 2-6) - є (Х, 2), где Аа) - ПОР) - 0 принадлежит заданному линеалу. В итоге 





42. (14) 


имеем 








вам рН" де Їж - но 
(е-Ю'я(ачН) ж 


22(8) ном с ОЛЯ 


Ле(х, с) з о гарно 


"Чароная АПВОЖЕХОТЬ 
Й З СВИСТ ПР САС ДЕ 
(а 





6 -хуз ОО НОГРО (ою не вДООГ У 
Отсюда 


го 29 -хуна (9 - ГО 
а 1 а 
Го во ЛОЮ їв б зарано я їз Ле г. зла 


2 Ге(е) я Уа (Та(5) чн в(БУНООУЇ - 35 - х)") 
ТСВРУЖИ МОВА ТО ВМ 


Алгоритмь. На основаниий (14-15) можно построить несколько способов восстановления 
контакта 2-2(3), асхдехЬ, по заданному полю 2(х), ссхха, |а,Б| с |с, |, 


зі пає. 














п'(2є. (15) 


пає8|а, Б) « тез | с, 4|. 
1. Способ Ньютона. Считая известньми а,Р, Н , разложим в ряд (14) относительно 2(С), 
з первьтм членом. УчитьІваЯ, что 
НІ Са хуя Н' | (б-Ю'яі  (Є-х) НО ру 22 


от (Е-хузо (Є-х) путом ГЕ-ху ча РТ з (вою 














находим 
Іа оман? 2«(а) 
з| да(54є. 
ко сон ово 
Отсюда, полагая Д 2(2) є 2,4 (2) - 2, (2), имеем 
г 22, (5) б -Н? 
; Ат СУ 45 - 8(х) як 
озарезе  сн ік ар 
я «орав (6) 





Ц -х)' я я, (5) 


2 
2 Способ "контакт 1". Равенство (15) запишем с точностью до Й (2), т.е. 


УДСО Б НИ 0 со рда С ЛИ, 
Гачок ою Ге ккОчнОЇ "У 
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2-х) -І2(Е)-О(КУНОЯЇ |з 
Є 2 232 П (Дає Б 
" 8-Х НГ е(а) (НОЯ) ) 
Сложим полученное тождество с (14), зано предварительно в (14) 2(х) в правую часть, 


б РО 
Гах Ма - їв поч НОЇ 



































З ЕКЗ ЕНН рою 
йо ПачхнаФНВФО кота Роби 
Тогда, обьединяя первьій и последний отвать стоящиє пе получим: 
24(5) -д'я Н? 
еріче Ф ние Й ЕР РОВУ ЛОТ В і; 
б: «Пе 8(ДНОУЇ з 
- 8397 КО ЕНСОМОТ о 
кон Па-юзнаДвткОЙ 
В зтом равенстве, как и раньше, б В итоге с точностью до Й (С) находим 
2:(2) -ху-Н? 
а а й 17 
Гаю або їв нина ко їі 
Полагая 2, (4) з 2, (2)--Н, (2), п я. ) 2,..., и А (2) - 0, приходим к уравнению 
ро 44 мо сад ро Крани з Я 
Гак на ис от у бо вбо) чи 


Оно разрешимо при известном с(5): 25) 7 2005) Йо(о) - Так как Й, (5) - що(5) 5 Зра(5)» 
то 





22, (9) (Є-ю'яН? 
нів з азоту - 
о. Ре рана еопеори 
1 2 
Гатне ніде аву 


откуда очевидно: алгоритмьт Ньютона и "контакт 1" формально совпадают только на первом 
шаге, так как 2, - 2, ана всех последующих шагах они формально различньі. 


3. Способ "контакт 2". Сложим величиньт А и 1. и, с точностью до НО , получим 














22(6) вва; 
4 
озеро: (5) ба Їю сс б МРТ аЄч 
ху? -2Д са їй (Є-х) ЕН? з 
Чо хг ОЇ (є) 2 (є-ху за (9 ачче(х). 


Отсюда, обьединяя первьій и последний интеграль! в правой части, а также полагая 


зма ФО, пед,1,2,...ив(Д2е0,ахсєхь, 
22, (5) 











приходим к уравнению Ме б кт паФЄ З 
оз - Н? бод ра 29 й 
і ас п й а й 19 
-Їн ти м заД 209)ч а ху ОР , (5)4с (19) 
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Зто уравнение, как и (18), можно решать итеративно, если задано 2, (2): 2 (4) 2 2606) П, (2). 


Уравнениє 0 (19) также | можно дать в виде,  аналогичном (18):  полагая 


п (2) з 2, (2) - 2, (2), получим 


22, (8) еоюонН" 








са М а ; а6- зе 
Їс -х)? зга 2:45) Си (5) ба їв -х) 8 2 (7) 2 200) 
22, (5) -х) й 2 - 
: 1 а 2 б 2 Гелі "би аб. 19 
Гео чні? сю Я а нд Ра 9 


А что, если определять не только контакт 2 - (2), ас 2 «Р,но и "асимптоту" Н з сопзії ? 
4. Способ Ньютона. Примем во внимание, что 
д вони 24 03 К8- хуяНн"' ОН 


ді (2-хуча (ре но з! за др» цеоЮзНнО 


и запишем 
(Е-хуя Н" 22(0) 2Н 
ЕІ ДЕННЕ В ДЕННЕ а зара ви у 
іє Йти у Мо б і 
Отсюда, полагая ча АЮ я НаяН,-АН,, пад, !, ій и А2(2)2-0, 
ДН,е0,ахєс с «ь,изадавая Н, и 2,(2), получаєм 
| 22, (5) 
по оо 











аг Дає о Н,4Є- 


5 ; Д 

ун; 

-Їи (Е-ху'яН; 
(Е-ху ча (2) 


Зто же уравнение в другой форме (непосредственно в "приближениях") имеет вид 


45 - 8(2. (20) 











го 22.Д со Он г. (2-Ю24Н? 
«ра ну НЕ а Пп руль 
Ган ногу Не о а а 


22, (5) 
(3-х за (2) 


5. Способ "контакт 17. В отличис от п. ра здесь помимо определения 27 ахісар, 
р п 








со Он 
-вооч| 2 Дає уні Н, 482, пе0,1,2,..., 00) 


будем находитьи Н п ФО Є «5. Для зтого рассмотрим функционал 
(а-х'-н? 


4 
созноЮ ; 


и определим на линеале его приращениє ДЛ 2(х, 2, Н) - є(х, 2-Й; Н --т)-- є(х, 2, Н), где 





2(х, 2, Н) «Й 


пса) « АС) «0. После несложньх за в итоге получим: 


22(5) 

суч фай Ї п 

й под'я(н жо)? Бо нан | 
он СРО 

зраннааткс - х) -(Нбт)") Раз 
ІЗб-юя(Невтуї 





гає 





а 








ов х'ян" р 
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2-х о ДФ з 2 (с 9 ПУЕ «(алу 
1 а 
«о ху Є | Е-'-н(ан Оп ГТ 


2 
Сложим полученное тождество с точностью до Й" (2) с равенством (14), перенесем 2(Х) в 








542. (21) 


правую часть полученного соотношения и сведем логарифмьг: 


сор нн 
Гаунтд ОЕЧГО 





ун? таб 














СГ (е-хуа(Нат? 0 (ох а)? г. (2-ху зн? 
І п а І а 
откуда 
М 2:(5) 2Н -хуз(Нату 
азу М заоч Че Їм | -. нон) 447800. 


Полагая в зтом уравнений 
зт (5) є з), На є Н, кт, пе0,1,2,..., и (220, се0,ахєхь, (22) 
и, задавая Н., 23), а 2 2 СБ, получаєм уравнениє 














в і рані (2) -х) ня Ні 
| с най Ма" са е їн РИ уодіб во- Оз 
Так как й, (2) 2 2,,2(4)-5, ми гами он Н,, » то уравнение (23) примет и такой вид 
р 2 у венЮвнноЗ 
2 2 Ной Ба І т а 
Гек над за іс і нн т а о 
р да 
ча у (6) уні ні 012» 23 
Час уні? б УНІ Нана п 5) 


6. Способ "контакт 2". Сложим величинь (14) и (21) с точностью до (2), и получим 




















М 2: (6) г В ою чіН ат)? й 
азу на Пе прнан 
уча 1 
40. 
ян нний 8" нене очної є 
Отсюда, полагая (22) и, задавая Н, - я 22200(0,а зу. «Ь,будем иметь 
г 9510) ЕО ен; 
і а со г п ен НЕ 
з аеваи З БАРАН ї зноЗнео Їте про 
5) 
42. 24 
о Пе З чні б ної пе 1 п а ча 
Подставляя йа(5) З но о) з (5). Тен - Н,о - На» получим то же уравнение в "при- 
ближениях" 
ро 28) ран 
5 2224 За Н,о ох б 
Газа! що . ху'яН; іш Ї» з уча СУ - 
да) 2Н Й 
- рот а па еооН 1 4б6- 
509 Ре з) все а 45 24) 
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2-х) Н: є-ху «ОД ; 
п -Н 4 п й - асе, 1 
«ча о - НІ) сб а , б є 7 Ї(2- за р сСаур Бо (5) ра (21 п 


Вьбор нулевого приближения. Пусть кривая 2(х), сс х «а имеет на |с, «| несколько 
зкстремумов в точках с, є |с, |, і - 1, 2,..., п. Так как 6.» 0, то каждому зкстремуму отве- 


част однойменньй зкстремум 2 - 2(2), а 2 4Ь. Определим во множестве (ехіг 2(х)) най- 
Іс, 4) 


больший и найменьший зкстремумь (наийбольший максимум с; и найменьший минимум 2,), 
т.е. 
23 З 8(с;) з 8ир ех 809 8; з Є(с,) з Бех 20). 
С, С, 


Обозначим через 2; - 2(с;) - 8ар ема с(6)), зо сс,)з їшбехіг с(С2)) соответствующиє 
(с, Ісаї 


зкстремумьі с. Здесь абсциссь 8 иа (и 8 2,) совпадают, что не всегда имеет место, но в 


данном случає верно. Нулевоє приближение контакта можно задавать, по меньшей мере, двумя 
способами (Ох направлена вверх) 


5 п 5()2 ЗНА, Но ахєхь, (25) 
ого50)2-НонкЛе(2),ахебєхь, (26) 
где Н,» 0 - задаваємая константа - "асимптота" функций 27, на левом конце промежутка 
Іа, БІ, Дде(2) «є(2)- 2(а), Де(5, Н,) - аналитически продолженноє значение Д2(2) на 
уровне - Но, Ко - КОН.) - подлежащий определению козффициент. 
Рассмотрим, прежде всего, способьт определения Ки. 
1. Пусть в точке с є |а, Р| располагаєтся 2, или 2,. Очевидно 
Ь 
зу Пи сан; Ре и и'-Н? о с 
22 (2-су с (О 5 ис пас 
где йикб-с, дфизаб, аза-с, Бебо-с. Возьмем 2, 2 2(2), (2, 4 2(2)), ачслЬ 
тогда, 





ші СЕ - су на Орел СЄ з с)" ча 1, (2 о (2012 (Є су чо І. || (27) 


Отсюда следует 











Й 2 2. Б 2 2 
и"-Н и"-Н 
8; 2 Пп но ди, ера г ди. (28) 
а 5 а і 
яз 2 
Оценим интеграл і гру В я полагая, что а-0, Б-і ар--с, Б -1-с. После 
2 іо 
ряда преобразований окончательно находим 
Ь 2 2 2 2 2 2 
и"ЯН а-су ян сЯанН 
Па-5у «2 (1-Ф п русі 
аа НО П(-су а с'Яа 
ра 
2Н агсі 22 агсі руч 
нео вк о 
Позтому 
1-с)-нН? 2-Н? 
РВЕ пн Лв Р и й 
19 па оо цк УА поро иа 
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-2Н ав ПЕ Р авто су (29) 
с Пс) об феої) 
«а п (- с) -Н? Р 
п с.) б. Ж 
8 ) Р а / сут; 
Н 2; 
ж2Н нон ДАНО (30) 


І 
Рос (1-с,) 2, - с, 1-с,) 
олученньюеє неравенства являются исходньютми для определения 2. ий сх. при заданньвах 3 
Пп 5 Н 
8; 171 8; . Отьшскание значений ра и ран однотипно, позтому в дальнейшем индексь опускаются. 


Составим функцию 


























тооза-сшіа- сна 1нсіс ча ага 
є 
Н 
-а-сушіа- с) - Н'|- сс '- Н')- а тузарестеноци 2. (31) 
-с(-с) 
Вьчисляем первую производную 
22 
«дай-ос с -Загсіє-з--- 
о Гатдічо сія ой -с(1- Й 
22 2-са-д-22 
т 22 2 2 
2 (с-с(1-с)) 
Петер п 
27 -с(1-с) 
2 
«(1-с) - РО ца 2, З дагсіїв-, Е ек ем 91 (32) 
(-с) ча са агезе (енд) Ка -єСреау ва 
Для вьічисления /"(2) производную /(2) лучше представить в виде 
х Ра 1-с с 
Кгу-ХІ- є -О2с -Загсіє ---- -агсіє - - 
ПоОех адіно сад Ше 9, 
геї 2а-да за) 
(2-са-Фуча (п'-с(-суча" 
1- 
поскольку агсіє ав -рагсі є не агсі є ПЕ РЕ . Дифференцируя второй раз (32) по 2, по- 
ї А 2 -с(1-с) 
лучим Же 
44-сг ст 232 с -с)) 
(а-с'я) (са) (г-сО-сФ) ча 
пе (22 -2с (1- суч (о с -с)) 33) 


«с зведено 
Приближения с, какого-либо из значений с. и 2, определяєм по схеме Чебьшева |Бере- 
зин, Жидков, 1960): 
са) РОС 
Син - Си ЛП І М у (34) 
ле . ХХ 


0 
Последовательнье приближения корня 2" уравнения (2) - 0 можно получить значите- 
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0 
льно проще, чем сделано П.Л.Чебьшевьм. А именно, если корень ГА ; располагаєтся в достато- 
чно узком интервале сі, 1 , ТО, ЄСЛИ 2, вьібрано с того же интервала, с точностью до мальях 3- 


го порядка относительно 20) на 19, 2, Є, ф| имеем: 
бо (0 д , (0) "Са, (0) 2 
з (с )2 (2), Ха чем ен 


Отеюда при /(2, ) я 0 (зто может бить, ибо только для кратньх корней (2) - Г'(2) - 0) по- 





лучаєм 
(0) з Ло, ЩО, (0) 2 
аобенк ен р Ск о 
Ре ЗРК) 
Учитьтвая малость разности (2) ні! РА ,ь на оснований предьщдущего соотношения положим 
2 
2 
: "ДРО о ІДУ ЕН ВЛ 

оо ЕЕ ся откуда Р позтому, 


РАС М Р З СИИе 9) 


г 21-Ї0) Раді.) 
"ОТ Ра) ЩО 


сходится (3-4 итерации для получения 20) с точностью до 10, что согласуется с результата- 
ми, полученньми в работе |Березин, Жидков, 19601). 





, т.е. получим (34). Процесс (34) при удачном вьборе 2, бьшстро 


Замечание І. Из (29) и (30) вьтекаст, что /(х.)20,а /(2)20. А при не очень мальжх 
Н » 0 имеем строгие неравенства уравнений 
Леї20, пе2«0. (35) 
Позтому корни 2,, 2, СОоответствующих уравнений /(2,, 2:)- 0, /(сць 2;)- 0 дают 
приближения 2. и 2, снизу и сверху соответственно, что может сказаться на "качестве" нуле- 
вого приближения 2, 7 2,(0): 
Замечание 2. Согласно замечанию І 2. и 2, определяются не очень хорошо при больших 
Н , особенно, если сами значения с, и 2, достаточно большие. Позтому дополнительно можно 
определить, если 68 - є(Р)- 2(а) є 0, кроме 2, и 2,, заданной левой Н' еще и правую 
Засимптоту" Н" кривой 2 2 2(2). Ото делаєтся по описанной вьшіе схеме, только в вьтраже- 


ний /(2) вместо є, или 2, берется б є для определения Н". 


Вьчислив так для каждого заданного значения Н,, і - 1, 2,..., п, левой "асимптотьт" ис- 
комой кривой 2-2(3), ах, определим значения параметров 2 . о о СН 
і-1,2,...,п, 

ид з20-НІ по та-НІ, АН, еНІ-НІ 
5 і? п) і і? і" і і? 


б 8 - 8 (Р) - є (4), б 8 ж 7 85 7 84), бЯ іа 7 8: 7 Є (а); 
и найдем 


у со 
пла ш ра б (36) 
ізі 








( 
Ф -і и , ЕФ Ро Ма 5 їз дні 
ТО о 3 би 7 Я 


Рассмотрим схему отьскания нулевого приближения для контактов, описьваємьіх функ- 








тах 
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циями из классов "антиклиналей" Оз (а, Б) и "синклиналей" бу (а, Б). Типичньми 
представителями зтих классов являются функции (рис. 4): 

а) антиклиналь 


«2)2-Нчай1і-с)8, 0«6х1, (37) 
6) синклиналь 
«2)2а-Н--ай0-2)6, 0«сєх«1. (38) 


Для отьскания К, систему координат удобно вьгбрать так, чтобьт точка зкстремума 2(х) 


находилась в начале координат, т.е. ехіг 2(х) - 2(0), тогда и интервал |с, | вьбираєтся 
Іс, 4) 


"симметрично" относительно начала, так же каки а--с, Рис. При зтом для модели (37) 
(слева на рис. 4) имеют место соотношения а, - агсіє С/Н , На, « 0, (2 8», а для моде- 


ли (38) - соотношения а, - агссє С/ 2, На, » с, а «5. 
Функция /(2) - 2с (соч) 4 агсів 3 -2с ш(с'-Н')-4Н асо где 8-8; 
1 


- в случає антиклинали, 8 - 8, - в случає синклинали (оба значения - со своими знаками). Не 


4с 


та зі . Так как Г(Н) - є, то в случае синклинали 


производнье /(2) - 4агсіє -; "есе 
2 С 


4 
ДОН) «0, ав случае антиклинали Г(Н)»0. В обоих случаях /"(Н)- по «0. Поз- 
Св 


тому Г(НУГ"'(Н)» 0 для синклиналии Г(Н)7"(Н) «0 для антиклинали. В обоих случаях 
для определения с, и 2, используєтся процесс (34). Попьштки обойтись грубьшми приближения- 
ми 2; й 2;, вьічисленньми не по (34), вообще говоря, успеха не имели", хотя в ряде случаєв 


приводили к хорошим результатам. 





Рис. 4. Модели контактов, описьваємьх функциями "антиклиналей" и "синклиналей". 
Рассмотрим некоторьте из попьмток. Пусть вначале 
2 2 
Сян с с 
8 -22 т - о НАН агсів -- -4г агсів-. 
СВО Н ра 
Для синклинали (рис. 4) с-Нчи, Н, 2, и» 0, а, 2. а, (знак равенства достигаєтся при 





с -» 00) и, следовательно, АН ав -4(Н чи) агсіє «ин 4На - «Н чи)а, «0. Позтому 
и 


" Слабая сходимость (до 50 итераций) и низкая точность, численная неустойчивость. 
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сЯН. є -5. ав 
а з са к(саН)е г, сс Н'уе 2-й -Н. (39) 
о 20 


Для антиклинали 2-Н-и, Н,и»д0, а, З Ц» и, следовательно, 


То 


4Н ав КН -и) агсіє - -4На -4са, 20. 
ус 


Позтому 





2 2 8 
сн ча 
ррайнальк А 3 МА РИН с (40) 
стЯа 2с 
При небольших значениях и (по сравнению с Н ) можно поступать по другому, а именно, 
принять 


2Н 
шіс' -(Н --и)/)-щ(с'-Н') я ТА 
Ге'нСН му Ну сут 


б 


с С с сн 
Н -ки) агсіє ----- Н агсіє - я» | агсіє -- - ------- |М 
ой б ЯН пт 


с 
для синклинали. Позтому 9 --4 агсіє Н и и, значит, 


з к/я аа (41) 


поскольку для синклинали 9 « (0). Для антиклинали, хотя 2.2 0,но 2- Н - и. Позтому фор- 


мула (41) остается прежней. Заметим, что Г(Н)- є, Г(Н) «4 агсіє - и позтому (41) явля- 


ется первьшим приближением в итерациях (34). 
Второй способ основан на определениий различньх моментов кривой 2 - 2(2). В просте- 


о 
йшем случає, по теореме К. Гаусса Ї 2 сах «2 т (т - возмущающая масса, сосредоточен- 
-о 
ная в области, ограниченной кривьми 2 з 2(2), 2 --Н, ас с, ось 02 направлена 
вверх), позтому 


та Ісав я | снає - анна, с-і. 


т 
-Н или 








т 
Отсюда 2; 2 Ні 
Ь-а ь-а 


1 о 
зро цу /8О0ахН, рр 


за РТ Ї еодах-н. (42) 


2л( 


а а 
В последних формулах полагаєм Ї (мах «| 2(хах. 


- і 


(Е-ху-Н? 


а а Ь 
Обратим внимание на  вьіЧИСЛенис 2( хх а | ах Пп 42. Подставляя 
| | | -ю'нФ) 


изб-х, дис -, ие ос, и, с - 4, получим 





ї б и'- НН? и и)" 
Ї олако Ди Р хана Зо аб- 
Я - иа Н ра 


и» 
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ь 2 2 яте 
(4-94 Н (5-9) (2-2) 

- іп 2Н агсів39--12-2 2231 42 
| еодто чад ЗНОВ о НОВ і З 
«сію сон аю ЕХ у2Намір о - 

й ся (2) азс) Н 


4 С 4 
2Н агсіє гла 2:(2) агсіє «9 -2:(С)агсіє Р ас- 
б Ї аа?" чн а"(2)1- сосна (21 баків б -2г(б)агоів 4 
б «(6) (5) 


-(ь-а) (сте зн?) -аша? я Ну ОН асо 2Н агсіє 3 


О численном решений. По предложенньм алгоритмам составленьт программьт и проведено 
численное моделирование. Теоретически зти конструкции показьтвают несколько лучшую 
сходимость при вьтчислений способом Ньютона по сравнению с классическими примерами 
ІГравиразведка, 1990), но их практическая зффективность требует проверки на подготовлен- 
ном материале полевьх гравиметрических сьемок. Попьтки обойтись грубьшми приближения- 
ми с, й 2, вьічисленньми не по (34), успеха не имели (крайне низкая сходимость, на поря- 


док меньше заявленной |Березин, Жидков, 1960|), хотя иногда приводили к математически 
корректньм результатам. Однако, геологическоє содержание грубьїх приближений едва ли 
достоверно. 
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Ю.Й. Дубовенко, О.А. Черная 


О вьюборе нулевого приближения при определениий сложного контакта 
Ю.І Дубовенко, О.А. Чорна 
Про вибір нульового наближення при визначенні складного контакту 


Ключові слова: теорія потенціалу, інтерпретація, аналітичні моделі, контактна задача, класи контактних 
поверхонь, розділення гравітаційних полів. 

Ключевьюе слова: теория потенциала, интерпретация, аналитическиє модели, контактная задача, классьт 
контактньїх поверхностей, разделение гравитационньтх полей. 

Кеууогав: рогепца! Реогу, іпіегргегацоп, апаїуціса! птодеіз, сопгасі ргобіег, сіа55е5 ої сопгасі 5игіасе5, єгамігу Пе!д5 
зерагацоп. 


Предложень аналитические модели гравитационного поля и геологической средь для горизонта- 
льно-слоистой средьт с несколькими плотностньми границами раздела слоев в нижнем полупростран- 
стве. Первая модель представляет собой уравнение силь тяжести, полученное из рассмотрения норма- 
льного потенциала силь тяжести в локальной окрестности текущей точки измерений в двумерном 
представлений. Для модели средь!, в дополнение к известньїм в теории интерпретации потенциальньх 
полей основньм классам контактньх поверхностей (Нумерова, Остромогильского и Страхова), введен 
новьій класс контактньїх поверхностей -- класс Черного. Приведеньт образць теорем разделения полей 
для случая нескольких односвязньх обьемов и нескольких непересекающихся слоев (два самьгх попу- 
лярньх начальньїх приближения средьт) в изучаємой модели. Описан новьій численньй подход к оп- 
ределению начального приближения плотностного контакта и его асимптот в многослойной геологи- 
ческой среде с помощью нескольких 'контактньх" алгоритмов. 
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Рис. 3. К разделению полей от нескольких слоистьтх источников. 
Рис. 4. Модели контактов, описьваємьх функциями "антиклиналей" и "синклиналей". 
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